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Abstract 

We consider the motion of a vortex in an asymptotically homogeneous condensate bounded by 
a solid wall where the wave function of the condensate vanishes. For a vortex parallel to the 
wall, the motion is essentially equivalent to that generated by an image vortex, but the depleted 
surface layer induces an effective shift in the position of the image compared to the case of a vortex 
pair in an otherwise uniform flow. Specifically, the velocity of the vortex can be approximated by 
U ~ (h/2m) (yo — V^C) \ where yo is the distance from the center of the vortex to the wall, £ is 
the healing length of the condensate and m is the mass of the boson. 

PACS numbers: 03.75.Lm, 05.45.-a, 67.40.Vs, 67.57.De 



1 



I. INTRODUCTION 



Vortex dynamics in a Bose-Einstein condensate (BEC) has been studied intensively, ini- 
tially in the context of superfluid helium and later in dilute trapped BECs. The motion 
of vortices in both uniform and inhomogeneous condensates has been the subject of many 
theoretical works, and extensive reviews of these efforts have been given in ^,2]. 

In this paper we consider the problem of vortex motion in an asymptotically homogeneous 
condensate in the presence of a solid wall where the wave function of the condensate vanishes. 
Recent discussions (see, for example, |3( and references therein) on the motion of vortices near 
the surface of trapped condensates have questioned the relevance of the method of images 
in describing this motion. In that case, the nonuniform condensate density is approximated 
by a linear function that vanishes at the Thomas-Fermi surface, and the vortex motion 
can be considered to arise principally from the local density gradient. Here, we consider a 
rather different situation, in which the condensate density approaches its bulk value within 
a healing length £, and the vortex is located in the asymptotically uniform region. In this 
latter case, the local gradient of the condensate density is very small. The motion can be 
interpreted as arising from an image, but the depleted surface layer induces an effective shift 
in the position of the image in comparison with the case of a uniform incompressible fluid. 

Our geometry is two dimensional, with the vortex aligned along the z axis, parallel to 
the surface of the wall. The dynamics of the time-dependent BEC in the presence of the 
solid wall at y = is described by the Gross-Pitaevskii (GP) equation 

-2i^ = V 2 ^ + (l-|^|>, (1) 

subject to the boundary conditions 

ip(x,y — 0,t) — 0, < y < oo, \x\ < oo, (2) 



in dimensionless units, such that the distance is measured in healing lengths £ = h/ \/2mgno, 
where g is a two-dimensional coupling constant with dimension of energy times area, m is 
the mass of the boson and n Q is the bulk number density per unit area. Time is measured 
in units of m£ 2 /Ti and energy in units of h 2 no/m. In our units, the speed of sound c in the 
bulk condensate is c = l/v2. 

In the absence of vortices, the exact solution of for the stationary state of the semi- 
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infinite condensate is 



/(</) = tanh(|//V2). 



(3) 



In classical inviscid fluid dynamics with constant mass density p, the relevant kinematic 
boundary condition at a solid wall with normal vector n is 



where u is the velocity of the fluid. The corresponding problem of a vortex moving parallel 
to the wall is solved by placing one or more image vortices in such a way that condition 
is identically satisfied. 

For the dynamics described by the GP equation (JTJ), the density p oc \ip\ 2 is no longer 
constant, but rather vanishes at the surface of the wall. Thus condition is automatically 
satisfied, and all components of u can in principle remain finite on the boundary. Therefore, 
it may seem that image vortices are irrelevant in the case of the GP equation, so that the 
vortex should remain stationary away from the boundary (where the fluid density is constant 
apart from exponentially small corrections). Our numerical simulations show that this is 
not true. In fact, the vortex moves parallel to the boundary, and it moves faster than a 
corresponding pair of vortices of opposite circulation in a uniform condensate in the absence 
of the depletion caused by the boundary. The purpose of our paper is to study this motion 
in detail. 

The paper is organized as follows. In Sec. II we find the family of disturbances moving 
with a constant velocity along the solid wall by numerically solving the Gross-Pitaevskii 
(GP) equation in the frame of reference moving with the disturbance. In Sec. Ill a time- 
dependent Lagrangian variational analysis is used to find the first two leading terms in the 
equation of the vortex motion in the limit of large distance from the wall. In Sec. IV an 
alternative approach based on the dependence of total energy and momentum on the vortex 
position is used to determine the vortex velocity. In Sec. V we summarize our main findings. 

II. Numerical solutions 



In what follows we seek solitary- wave solutions of Eq. (0) that preserve their form as they 
move parallel to the wall with fixed velocity U. For each value of the velocity U, we have 
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p u • n = 0, 



(4) 



FIG. 1: [Color online] Graphs of the velocity of the vortex U versus the vortex position yo a s calculated 
via numerical integration of (J5J subject to the boundary conditions tp — > 1 as x 2 + y 2 — > oo without the wall 
(solid line) and the asymptotics given by U = (2y )~ 1 (dashed line). 
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where r\ — x — Ut. The GP equation (|T} becomes 

2i£^ = V 2 V + (1 - M V, (5) 



where we set x = 77. In the a 
found by Jones and Roberts 
and energy E, given by 



)sence of the wall, the solitary-wave solutions of Eq. (jSJ) were 
For each value of U, there is a well-defined momentum p 



41 



p = k I W ~ l ) d *^ - (V> - 1)^1 ^ , (6) 

E =\ \ \Vi)\ 2 dxdy + \ J (1 - \i)\ 2 ) 2 dxdy . (7) 

In a momentum-energy pE plot, the family of such solitary- wave solutions consists of a single 
branch that terminates at p = and E — as U — > c (we call this curve the "JR dispersion 
curve"). For small {7 and large p and the solutions are asymptotic to pairs of vortices 
of opposite circulation. As p and E decrease from infinity, the solutions begin to lose their 
similarity to vortex pairs. Eventually, for a velocity U ~ 0.43 (momentum po w 7.7 and 
energy E m 13.0) they lose their vorticity (i/j loses its zero), and thereafter the solutions 
may better be described as "rarefaction waves" that can be thought of as finite amplitude 
sound waves. The velocity of the vortex pair in the absence of the boundary is plotted as 
a function of the position of the vortices ±t/ shown in Fig. The dashed line gives the 
asymptotic velocity valid for large y as U = (2y ) _1 - 

In analogy with these results, we used numerical methods to find the complete family of 
solitary- wave solutions of (0) subject to the hard- wall boundary condition (J2J). Specifically, 



FIG. 2: [Color online] The momentum-energy curve of the solitary wave solutions of Eq. (J5J with the solid 
wall (solid line) and the JR dispersion curve (dashed line) that has no solid wall. For the solid-wall boundary 
condition, the vortex solutions with nonzero vorticity and nodes are shown in black and the vorticity-free 
solutions in grey (green). 
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we mapped the semi-infinite domain onto the box (— f , f) x (0, |-) using the transformation 
x = ta,n~ l (Dx) and y = t&n~ 1 (Dy), where D ~ 0.4 — 1.5. The transformed equations were 
expressed in a second-order finite-difference form using 200 2 grid points, and the resulting 
nonlinear equations were solved by the Newton-Raphson iteration procedure, using a banded 
matrix linear solver based on the bi-conjugate gradient stabilised iterative method with 
preconditioning. Similar to we are interested in finding the dispersion curve for our 
solutions in the pE plane. The energy and impulse of each solitary wave are defined by the 
expressions ©-© appropriately modified for the "ground state" given by f(y): 

P= 3 / K^* - f(v)) 9 ^ - (V- - f(y))d x r\ dxdy , (8) 
E = i JOWl 2 +i(l-H 2 ) 2 -sech%/ v / 2)] dxdy. (9) 

In Fig. El we show the resulting solutions in the pE plane. The plot of the velocity depen- 
dence on the vortex position is given below in Fig. |3] in Sec. IV. All our vortex solutions 
with a rigid wall (those with a node in the fluid's interior) move with velocities less than 
U = 0.47. For U > 0.47, the zero of the wave function occurs on the wall only, and the 
solitary waves resemble rarefaction pulses of the JR dispersion curve away from the wall. 
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III. Variational approach 



The time- dependent variational Lagrangian method offers a convenient analytical ap- 
proach to estimate the vortex velocity for large yo- The dimensionless GP equation is the 
Euler-Lagrange equation for the time-dependent Lagrangian functional 

c = r - e = |i/ (r^t - dxd y - § / (iwi 2 + \W) dxdy (10) 

where the time-dependent terms constititute the "kinetic energy" T and the remaining terms 
are the GP energy functional £. 

We assume a trial function that depends on one or more parameters, and use this trial 
function to evaluate the Lagrangian L in Eq. which will depend on the parameters and 
their first time derivatives. The resulting Euler-Lagrange equations determine the dynamical 
evolution of the parameters. For the present problem of a vortex moving parallel to a rigid 
boundary with the boundary condition (J2J), the vortex coordinates (xcb2/o) serve as the 
appropriate parameters, where — oo < xo < oo and < yo < oo. 

When the condensate contains a vortex at a distance yo from the boundary, the original 
condensate wave function (|3J) acquires both a phase S(r, r ) and a modulation near the 
center of the vortex, where the density vanishes. To model this behavior for the half space, 
it is preferable to include an image vortex at the image position (x , —yo). In this case, the 
approximate variational contribution to the phase is 

S(r, To) = arctan ( — — — ) — arctan ( - ^° ] , (11) 

\X-XoJ \X-XoJ 

where the second term reflects the negative image vortex. The derivation of the Gross- 
Pitaevskii equation involves an integration by parts of the kinetic energy density |V-^| 2 
to yield —ip*'V 2 ip plus a surface term proportional to ip*d y ip and this is one rationale for 
including the image. Strictly speaking, this term vanishes because /(0) does so, but the 
image vortex ensures that the contribution vanishes even in the case of a uniform condensate. 
The image vortex also cuts off the long-range tail of the velocity, giving a convergent kinetic 
energy even for a semi-infinite condensate. It thus seems more physical to include the image 
in this particular geometry, even though the image is often omitted for the highly_nonuniform 
density obtained in the Thomas-Fermi limit for a trapped condensate 

In addition, the vortex affects the density near its core, which is modeled by a factor 
v (|r — ro|), where v(r) vanishes linearly for small r = x 2 + y 1 and v(r) — *■ 1 for r S> 1 P, ^. 



tie highlyn 
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In principle, the function v(r) can be taken as the exact radial solution of the Gross-Pitaevskii 
equation in an unbounded condensate, but this choice requires numerical analysis, and it is 
often preferable to use a variational approximation. A particularly simple choice is [l| 

{r/A for r < A; 
(12) 
1 for r > A, 

where A is the effective vortex core size; a variational analysis yields the optimal value 
\ = VE. With these various approximations, the variational trial function is fl 

V(r,r ,t)= e iS ^f(y)v(\(r-T \). (13) 

The time-dependent part of the functional in Eq. (fTU|) becomes 

r = -J [f(y)] 2 d t S \v(\r-r \)\ 2 dxdy^- J [f(y)] 2 d t Sdxdy, (14) 

where the last approximation omits the effect of the vortex on the density, replacing | v | 2 by 
1 throughout the condensate. A straightforward analysis then yields 

ryo 

T^27tx [f(y)] 2 dy, (15) 
Jo 

where Xq is the velocity U of the vortex parallel to the wall. 

The contribution to T from the vortex core yields a term that is smaller than Eq. (|15|) by 
a factor of relative order y^ 2 , which is negligible relative to the leading correction of order 
y^ 1 that we retain here. 

Since the energy functional will turn out to depend only on the single coordinate yo, the 
Euler-Lagrange equation for xo implies that yo remains constant (as expected from energy 
considerations). In contrast, the equation for y reduces to 

d dC dC 

dt dy dy 

since yo does not appear in T (and hence in C). Thus the dynamical motion of the vortex 
is given by 

2n[f(y )] 2 dyo 

It is evident that only the derivative dS/dyo is relevant, so that several terms in S 
play no role in the present analysis. For example, the derivative of the interaction energy 
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£int(z/o) = i / dxdy vanishes exponentially for y 3> 1 and does not affect the dynamics 
of the vortex for large jjq. Similarly, the kinetic energy in Eq. (fTUj) separates into two parts, 
arising from the density variation and the flow energy, respectively; the contribution from 
the density variation is also negligible for yo ^> 1. 

The remaining (dominant) kinetic energy Sk v — hi \^S\ 2 \ip\ 2 dxdy, arises from the vor- 
tex flow. The squared velocity now follows from Eq. (jllj) 

1 1 



\VS\ 2 = ^ 

x ) 2 + (y-y ) 2 (x - x ) 2 + (y + y ) 2 _ 
The resulting flow-induced kinetic energy is 

= f(y) 2 v(\v-r \) 2 \VS\ 2 dxdy. 

It is convenient to divide this integral up into three (strip-shaped) regions 

region I: < y < y — A (note that v = 1 in I) 
region II : y - A < y < y + X 

region III : y + A < y < oo (note that v = 1 in III) 



(18) 



(19) 



(20) 
(21) 
(22) 



In region II, inside the vortex core |r — r | < A, the integrals can be found approximately in 
cylindrical coordinates and yield £^ vc ~ (n / 2) f (y ) 2 , neglecting terms of order A 2 /y 2 ,. The 
remaining region of the strip II outside the core simplifies because v = 1. It is convenient to 
parametrize y = yo + Asin# by an angle 9 that runs from —tt/2 to tt/2. In this region, the 
symmetry in x allows us to consider only x > 0, and the lower limit for x is x m (6) = A cos 9. 
The relevant integral is 

1 1 



x m (6) 





1 


\y 


-yo\ 
\e\ 




sin^ 



_x 2 + (y - y y 



7T 

arctan 

2 



x 2 + (y + y Q ) 2 _ 

Xm(9) 



dx 



\y-yo\ 



y + yo 



7T 

arctan 

2 



x m (9) 

y + yo 



71 

4j/o" 

The total answer for region II is 



(23) 



kvll 



+ ln2 



A 

Wo 



(24) 



where In 2 arises from the definite integral $_^i 2 1^1 /I tan^| d9 = 7rln2. 
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In regions I and III, the integrals can be found in cartesian coordinates, integrat- 
ing over x first. Each of these gives two contributions; one is simply a combination 
of logarithms obtained by replacing f 2 by 1 and the other from the remainder with 
-(1 - f) = -sech 2 (y/V2). 



2 V A / 2 Jo \yo-» yo + y; 



21nl^Vln^ + A 



* tom " 2 r~ V a ; ~v a 

-J / sech 2 (y/V2) ( — ^- + —J dy. (26) 

To evaluate d£kvi/dyo and dSkvin/dyo, we first differentiate the expressions in Eqs. (J2*H)) and 
(J2HJ), expand the integrands in the powers of l/y through 0(yQ 2 ) and integrate to get 

dE ™ ~ 71 +^tanhf^V (27) 

(28) 



dy 2y -A y 2 \ y/2 

dSkvm ^ 7ryo 

dy ^ (A + y )(A + 2y )" 

A combination of these contributions gives the vortex velocity in Eq. (jl7j) as 



£7 = x ~ -coth 2 -^= 
2 v/2 



1 yo 



2j/o-A (j/ + A)(2j/ + A) 

+ v / 2^--^ + ln2Vech 2 (^)tanh(^) 
V 2 2j/ / K y/2 J y V2 ; 

+ — -^tanh 2 ( -^p=] + -^-tanh ^ "" 



2y 2 \V2~J yl \ V2 

This further simplifies to 

1 / 1 J/o A + 2y^ \ 



(29) 



C/ ^2l v 22/o-A + (i/o + A)(2 2/ o + A) + 2y 2 J' (30) 

after we approximate the hyperbolic functions by their large y$ behavior. When we neglect 
terms of order the expression (|3Uj) finally reduces to 

U*±(l + ^), (31) 



2j/o V Vo 

independent of A. 
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IV. Vortex velocity through the Hamiltonian relationship between energy and im- 
pulse 



In this section we present a different approach to the asymptotics for the vortex velocity 
based on the relationship between energy and momentum of the vortex pair. We compare 
the motion of a pair of vortices of opposite circulation (JR solutions) that satisfy 

2iU^ = V 2 ^i + (l-|^i| 2 )^2, W->1 M-xx), (32) 

with the motion of a vortex next to the solid wall 

2iU 2 ^ 2x = V 2 ^2 + (1 - |^2| 2 )^2, \*h\ -> | tanh(y/V2)| |x| -> oo. (33) 

For the asymptotics, we are interested in the solutions for small Ui, for i = 1,2, that 
correspond to a pair of vortices of opposite circulation. We calculate the following quantities: 
the position of the pair (0, =h/o)> the energy and impulse given by ©-© for i = 1 and by 
©-© for i = 2, so that 

Ui = -rr^. (34) 



rffn 

in |4|, |10J ; simi 



These expressions for 2 = 1 were derived 
the expressions for i — 2. 

Note that C/j, and Pj are functions of y and if y 3> 1 



similar arguments immediately lead to 



£7i = 27rlog(2j/ ), Pi = ^y , (35) 



(see, for instance, (llj). From (jHljl we have 

dpx/dyo 2y 



U. = = ± (36) 



as expected. 

For large y an accurate approximation to the solution of (32) for the uniform flow was 
found [l^ as ipi — u i( x ,y) + i> v i( x ,y) where 



Uifoy) = (x 2 + y 2 - y$R(y/x 2 + (v ~ y ) 2 )R(V^ 2 +¥+ 2/o) 2 ) 



Vifoy) = -2xy 0J R( v 7 ^ 2 + 5 - Vo) 2 )R( V* 2 + (v + 2/o) 2 ), (37) 

where fl(r) = ^(0.3437 + 0.0286r 2 )/(l + 0.3333r 2 + 0.0286r 4 ). Another accurate choice is 
R(r) = (r 2 + 2) -1 / 2 . Similarly, we expect that ip 2 is accurately approximated by ip2 = 
^!|taoh(j//v^)|. 
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The question we pose is: What is the position of the vortex (0, yo) moving parallel to the 
solid wall with the same velocity as the vortex pair at (0, yo — I) in the uniform flow? Thus 
we seek the solution of 

Ux(y - l{yo)) = U 2 (yo) = f^r^, (38) 

dp2/dy 

where we explicitly indicate the dependence of JJ\ and U 2 on the vortex position. Since 
Ui(y — I) = (2(?/o — 0) -1 j we obtain the expression for the shift in the vortex position, I, in 
the presence of the wall as 

,/ s 1 dp 2 /dy 

l{yo)=y °-2dE^dy- (39) 
We rearrange the right-hand side of (fHHJ) and use to obtain the final equation that 

determines l(yo): 

, , . 1 Arc + dp , s 

Km) = yo- (40) 

2 2Tt/y + dE 

where dE = d(E 2 — Ei)/dyo and dp = d{j>2 — Pi)/dyo in the integral form given by ©-© 
for Ei and pi and ©-© for E 2 and p 2 . In evaluating the contribution E 2 — Ei only the 
kinetic terms involving derivatives with respect to x were kept. The integrals dE and dp 
are exactly integrable in x with the use of Mathematica, in which the leading order terms in 
l/yo are given by 

g^. POO 

dp = 3 / sech 2 (y/V2)dy + O(y r °), 

dE = (- 2 + n{n2 ~ 18) ) sech 2 (y/V2) dy + 0(y^). (41) 

Wo z yo J j -00 

With f™ sech 2 (y/V2) dy = 2^2 we finally arrive at 

The vortex next to the wall moves with the velocity 

U 2 = ^—7^, (43) 

2(2/o-V2)' 

which is the main result of our asymptotics. Note that if we expand (|4*3*|) in a Taylor series 
we get U 2 = (2y )~ 1 (l + \/2/y ), which agrees with the result of Section III. Fig. El gives 
the plot of the vortex velocity U as a function of the distance of the vortex from the wall yo 
for the numerical solutions found in Sec. II, asymptotics found in Sec. Ill [see Eq. fj3~T)) ] and 
asymptotics (J4*3~|) . 
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FIG. 3: [Color online] Graphs of the vortex velocity U versus its distance from the wall yo as calculated 
via numerical integration of (|33|l (black solid line) and the asymptotics given by i|31fl (red short-dashed line) 
and by (|43|) (green solid line). Also shown is the velocity of the vortex calculated by numerically integrating 
the right-hand side of i)38fl (blue long-dashed line). As one can see the simplifications made to derive (|43[) 
are consistent with the full expression 1|38JI for yo > 4. 
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V. Discussion and Conclusions 



In a uniform superfluid with a solid boundary, the motion of a quantized vortex arises 
from the image that enforces the condition of zero normal flow at the wall. In a trapped 
condensate, however, the image is generally omitted. Instead, the motion can be considered 
to arise from the gradient of the trap potential, which is the same as the gradient in the 
density in the Thomas- Fermi limit 1]. This behavior is especially clear for a single vortex 
at a distance r from the center of a cylindrical container of radius R |6|. For a classical 
incompressible fluid, the vortex precesses at an angular velocity 

* =-m^—2 (44) 
ci m R A — 7"o 

because of the image vortex at R 2 /r . In contrast, the precession rate in a trapped cylindrical 
condensate in the Thomas-Fermi limit 

* HR/O (45) 



tf m R 2 — r 2 



is larger than (J44j) because of the (typically large) logarithmic factor. Although the denomi- 
nators of f)44j) and (|45j) both vary quadratically with ro, the first result arises from the image 
and the second from the parabolic trap potential (and thus the parabolic density). If an 
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image were included in the analysis of the trap, it would add a correction of order 1 to the 
large logarithm ln(i2/£); such a term is comparable to other terms that are usually omitted. 

As an intermediate situation between these two extremes, the present paper has analysed 
the dynamics of a vortex in a half space bounded by a solid wall on which the density of 
condensate vanishes. This geometry represents the simplest problem of a vortex in a conden- 
sate interacting with a surface. Since the gradient of the density vanishes exponentially for 
yo ^> £, only the image remains to drive the motion in the asymptotic region. Our geometry 
allows us to separate the effect of the surface from the effect of the density gradient, both of 
which appear in the more complicated problem of an inhomogeneous trapped condensate j^] . 
We found the complete family of solitary-wave solutions moving with subcritical velocities 
parallel to the wall. In addition, both a variational analysis and the Hamiltonian relationship 
between energy and momentum were used to give the velocity of the vortex as a function of 
its distance from the wall. These results are identical through to the first correction term, 
where the small parameter is the inverse distance from the wall. Our main results are (i) 
that the vortex moves as if there was an image vortex on the other side of the wall, which 
essentially replaces the boundary condition (J3J) with a more stringent requirement u • n = 
and (ii) that the depleted surface layer induces an effective shift in the position of the image 
in comparison with the case of the uniform flow. Specifically, the velocity of the vortex can 
be approximated by 

U « = , (46) 

2m(y -Vk) 

where yo is the distance from the center of the vortex to the wall, £ is the healing length of 
the condensate and m is the mass of the boson. 
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